/VJZ> 


The Astrophysical Journal, 488 : 881-894, 1997 October 20 

© 1997. The American Astronomical Society. All rights reserved. Printed in U S. A 


/W- 
c _ : 








tffc 7 




207118 

COMPTON SCATTERING BY STATIC AND MOVING MEDIA. I. THE TRANSFER 

EQUATION AND ITS MOMENTS 

p . . , . Dimitrios Psaltis and Frederick K. Lamb 

PhyS1CS aDd AStr ° n0my Depar ' memS > UniVe nl ,IHn °; S at L'rbana Chatnpatgn, 1 1 .0 West Green Street, Urbana, IL 61801-3080- 

demetns(aj astro. uiuc.edu, f-lamb(a) uiuc.edu ’ 

Received 1996 April 19; accepted 1997 June 5 

ABSTRACT 

Common scattering of photons by nonrelativistic particles is thought to play an important role in 

nh^tn n8 lr the t radiatl0n SP . eCtrUm ° f many astro Physical systems. Here we derive the time-dependent 
tbn L kme iC . equ ^ tlon that describes spontaneous and induced Compton scattering as well as absorp 
t.on and emiss.on by static and moving media, the corresponding radiative transfer eqL^on and their' 

me" um a We sho T "t m ° mentS ’ b0,h ta the system frame a " d the tome comovmg wl the 
order 7ht h ? / !?• neccssar >' to use the correct relativistic differential scattering cross 8 section in 

° h( f V t b a ph 2 ton T netlC equatl0n that is correct to first order in e/m e , T jm and V where e is 
the photon energy T and m e are the electron temperature and rest mass and V is the elect ron bulk 

tion l^t 1 " UmtS ° f S* S r d 0f , light - We als ° demonstrate that the terms in tEe radiative transfer eSU 
fr - ha *i ^ second order m V sh ould usually be retained, because if the radiation energy densitv is 
““ y ar ® e ’ c ™ par f d radiation flux, the effects of bulk ComptSSS described by the 

erms that are second order in V can be as important as the effects described by the terms that are^rst 

filJnd F ’ CVe / n W o n V 1S Sma "' ThC system ' and fluid-frame equations that we derive are correct to 
When F Cr m i^ e ' ? Ur sy ^ em ^ rame equations, which are correct to second order in V may be used 
when V is not too large. Our fluid-frame equations, which are exact in V, may be used whin V -* 1 
Both sets of equations are valid for systems of arbitrary optical depth and can therefore be used in both 
the free-streaming and diffusion regimes. We demonstrate that Comptonization by The ^ electron bull 

motion occurs whether or not the radiation field is isotropic or the bulk flow converges and hEt it k 
more important than thermal Comptonization if V 2 > 3 TJm e . converges and that it is 

Subject headings: plasmas — radiation mechanisms: thermal — radiative transfer — scattering 




I. 1JN 1 KODUCTION 


active galactic ntfdei (see mZZ °‘ ^ m51 

isgsaii wmmm 

r d “ ,s ™ s has mo,,vaied -- “ uih ” s *-£ i**- re"ge‘„s 

uSbSST JSSSS Isat'umb 01 T ZCStr 0 S ,he ^ (see. e.g„ 

moments for static and ^“vinf m’edf^ and found r Z } T aVC reder,ved , the radia “ve transfer equation and its 

kinetic or radiative transfer equations, as we explain in § 2 TteS^X? £a7 “‘faneSLI^ Z 
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differential scauering cross section in order ,o obtain a bulk »doa“n 

In^n4'^ 

£S.T£ 'lions Sc'ond'order terms 

Yin & Miller 1995). If the terms that are second order in V are instead neglected, signmca.. 

photon kinetic equation and its moments. introduce our notation In 8 4 we derive the time-dependent 

in 8 3 we state our assumptions and approx, mat, on ylgte" moving media, the correspond, ng 

photon kinetic equation that describes spontaneou moments, in fhelystem frame and in the frame comoving with 

radiative transfer equation, and eir z kinetic and transfer equations, because, although it is usually 

Kair^OT 

«■>*« “f i"** *“**• in 5 5 we ” mman “ our 

transfer equation that descT^s absorptr c d transfer J on witho ut any corresponding absorption term (see, e^g 
Bbndford'&^ayne'lffXl^ifis'fundanientally inconsistent with thermodynamics and leads to a radiative transfer equat.on that 
has a different mathematical character Iron, that of the ,he,modynam,cally cons, stent equat.on. 

2. MOTIVATION 

The radiative iransler equation that describes scattering of ^ to^it 

on,, derivatives with respect to the spatial the angular 

this equation is nonlocal in photon energy a p velocity distribution of the particles, and the differential 

dependence of the specific intensity of the radia 10 j , t „ ain better physical insight, and facilitate comparison 

scattering cross section. In order equft.on into a partial differential 

moments of the specific intensity of the radiation field. calculate the spectra of X-ray and soft y-ray 

PiilllS^ppias 

Xro^qr LTm- be ^ and some subde points that 

must be taken into account if the diffusion approximation is used. 

2. 1 . Approximate Scattering Cross Section 

In order to obtain a transfer equation that is consistently accurate to find order m e/m ,, it “ 

To first order in e/m,,, the Klein-Nishina cross section is 

(1) 


do 

o 


TS [l + (, « 


■ /o) 2 ] i - 2 — (i - i 0 ■ i'o) 

L m e J 


where o, is the total cross section for Thomson scattering and f„ and f„ are the direction vectors of the photon in the electron 
ofr«S™rSSCt sSSS. anMtron velocty, introducng errors in the radiat.ve 
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transfer equation of order v/m e , TJm e , and V. Finally, use in the system frame of the Thomson or Klein-Nishina expressions 
for the differential scattering cross section (Pomraning 1973; Payne 1980; Madej 1989, 1991) introduces errors in the transfer 
equation of order V and TJm e , because these expressions are valid only in the electron rest frame. 

2.2. Importance of the Terms of Order V 2 

It has long been recognized that if the divergence of the electron bulk velocity is nonzero, the electron bulk motion causes a 
secular change in the energy of the photons (see, e.g., Chan & Jones 1975; Blandford & Payne 1981a; Fukue et al. 1985). The 
photon kinetic equations derived in these works do not include any terms of second order in V and predict that photons are 
systematically upscattered in energy by the electron bulk motion if the flow is converging. Flowever, as we show with the 
examples that follow, photons are systematically upscattered by the electron bulk motion even if the terms that are first order 
in V have, on average, no effect. Indeed, in some circumstances photons are systematically upscattered by the bulk motion 
even if the flow is diverging . 

To see this, consider a situation in which the bulk motion of the electrons can be described as isotropic turbulence in 
which the velocity correlation length is much smaller than the photon mean free path and the bulk velocity satisfies V • V — 0. 
If all the terms that are first order in V are included correctly in the photon kinetic equation, but terms that are second order 
in V are not included (Fukue et al. 1985), the zeroth and first moments of the kinetic equation will indicate that the electron 
bulk motion has no effect on the photon energy distribution on timescales longer than the velocity correlation time. In reality, 
however, the effect of the electron turbulent motions under the stated conditions is completely analogous to the effect of 
electron thermal motions, that is, the turbulent motions cause the mean energy of the photons to increase steadily with time at 
a rate that is proportional to the second moment of the turbulent velocity field. Hence, the electron bulk motion causes the 
mean photon energy to increase with time, even though the flow is not converging; in fact, in this case no terms that are first 
order in V have any effect on the radiation field, on average. If the turbulent velocity is high enough, the terms that are second 
order in V will cause the mean photon energy to increase even in a diverging flow. 

Origin of terms that are second order in V. — The origin of the terms in the photon kinetic equation that are second order in 
V can be understood by considering a cold (T e = 0) electron fluid with a uniform bulk velocity V in the system frame. For 
simplicity, let us analyze the scattering of photons in the Thomson limit in the electron rest frame. In this limit, the angular 
distribution of the scattered photons is backward-forward symmetric in the electron rest frame, so the average photon energy 
in the system frame after scattering is y 2 (\ — • V)) times the energy before the scattering (see Rybicki & Lightman 1979, p. 

198), where y = (1 — V 2 )~ 1/2 and the average is over the initial photon wavevectors k in the system frame. Consider first the 
case of a photon distribution that is isotropic in the system frame. In this case the average photon energy in the system frame 
after scattering is y 2 times the energy before scattering, so in each scattering the average energy of a photon is increased by an 
amount of order V 2 and hence there must be terms of order V 2 in the kinetic equation that describe this effect. In the more 
general case of a photon distribution that is anisotropic in the system frame, these terms of order V 2 almost always produce a 
secular increase in the average energy of the photons, but there are also terms of order V in the kinetic equation that may 
cause the average energy of the photons to increase or decrease, depending on the velocity field and on the angular 
distribution of the photons (there are generally terms of order V even if V * V = 0; see Fukue et al. 1985). 

Relative sizes of terms of different order in V . — The order in V of a given term in the various moment equations does not by 
itself determine whether the term should be retained in solving a given transport problem. This is because the terms of 
different order in V involve different angular moments of the radiation field, so the relative sizes of the moments must also be 
considered in determining which terms should be retained (see also Yin & Miller 1995). 

For example, in the zeroth moment of the radiative transfer equation, which is a scalar equation for the energy density of 
the radiation field, the terms that are first order in V involve the scalar product of the vector quantity V with the vector 
quantity defined by the radiation field, i.e., the radiation flux H (see eq. [34]). In contrast, the terms that are second order in V 
do not involve H but instead involve the scalar quantity defined by the radiation field, i.e., the radiation energy density J, and 
the radiation stress-energy tensor K lj (again see eq. [34]). In general, V i VjK lj is of the same size as V 2 J. Hence, if J is 
sufficiently large compared to H , the terms in the zeroth moment of the radiative transfer equation that are second order in V 
can be as large as the terms that are first order in V. 

In the first moment of the radiative transfer equation, which is a vector equation for the radiation flux, the situation is 
reversed in that the terms that are first order in V involve J and K lj (see eq. [40]), whereas the terms that are second order in V 
involve H and Q ljk . In general, V } V k Q ljk is of the same size as V( V • H). Hence, one might be tempted to neglect the terms of 
order V 2 in the first moment of the radiative transfer equation if J is large compared to //, even if the terms of order V 2 are 
retained in the zeroth moment equation. However, this is in general unsafe, because it involves treating differently the same 
terms in the radiative transfer equation, from which both moment equations are derived, in computing the zeroth and first 
moments. In addressing a given transport problem, one can only determine which terms in the expansion in powers of V must 
be kept by considering the boundary conditions as well as the transfer equation. 

The radiation field quantities that appear in terms involving only odd powers of V are all of about the same size. Similarly, 
the radiation field quantities that appear in terms that involve only even powers of V are all of about the same size, although 
they may be much larger or smaller than the radiation field quantities that appear in the terms involving odd powers of V. 
Hence, when V is much less than 1, terms of order V 3 and higher may be safely neglected in the derivation of the transfer 
equation (see also Yin & Miller 1995). 


2.3. Use of the Diffusion Approximation 

In situations in which the change during a scattering of the energy of a photon as measured in the fluid frame can be 
neglected and the longest photon mean free path 2 max is much smaller than the smallest length scale L min on which physical 
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variables change, the specific intensity in the frame comoving with the electron fluid If(e f ) is given to lowest order in A max /L min 
by (see Mihalas & Mihalas 1984, p. 457) 

If( e f) = f * , ( 2 ) 

where J f (e f ) and H f (e f ) are the zeroth and first angular moments of I f (e f ). The diffusion approximation (sometimes called the 
zeroth-order diffusion approximation; see Mihalas & Mihalas 1984, § 97) consists in assuming that lj(e f ) is given exactly by 
equation (2). Then the first fluid-frame Eddington factor f l /(e f ) = K l J(e f )/J f (e f ) is exactly equal to ^ S lJ ; here K l }(e f ) is 
the second moment of l f (e f ). As is evident from equation (2), if the source function involves only the zeroth, first, and second 
moments of l f (e f \ then in the diffusion approximation it is only necessary to solve the three equations consisting of the zeroth 
and first moments of the fluid-frame radiative transfer equation and the closure relation f l f(e f ) = |c> tJ for the three moments 
H f (e f ), and K l j(e f ). The specific intensity may then be calculated in this approximation using equation (2). 

The moments of the specific intensity in the system frame do not satisfy equation (2), even when A max /L min <| 1. If equation 
(2) is mistakenly used to relate the specific intensity to its zeroth and first moments in the system frame (Blandford & Payne 
1981a) rather than in the fluid frame, then errors of order V will be introduced in the radiative transfer equation and its 
moments, even if A max /L min < 1 (see Fukue et al. 1985). 

If one wants to use the diffusion approximation (2) to solve for the moments of the specific intensity in the system frame, one 
must solve the four equations consisting of the zeroth and first moments of the system-frame radiative transfer equation and 
the two closure relations involving the first two Eddington factors for the first four moments J(e), K lj (e), and Q ljk (e) of the 

specific intensity 1(e) in the system frame. The reason is that when equation (2) is boosted into the system frame and the terms 
that are second order in V are retained, as in general they must be (see above), the third moment of the specific intensity Q ljk (e ) 
is introduced into the closure relation. (Fukue et al. 1985 were able to set up a closed system of equations consisting of the 
zeroth and first moments of the system-frame radiative transfer equation and a closure relation involving the first Eddington 
factor only because they neglected all the terms in the radiative transfer equation that are second order and higher in V.) 

When, as in Compton scattering, the energy of a photon in the fluid frame changes in a scattering, equation (2) may not be 
accurate for all photon energies, even if the mean free path of a photon is independent of its energy and A max /L min 1. As an 
example, consider a slab that is infinite along the x- and y-axes but finite along the z-axis, in which the electrons are cold and 
have no bulk motion. Suppose the slab is illuminated from one side with monochromatic photons of energy e in propagating in 
the z-direction. Since the electrons are static, any photons that have scattered have lost energy to the electrons and hence no 
longer have energy e in . Therefore all photons anywhere in the slab with energies equal to e in have never been scattered and are 
still propagating in the z-direction. As a result, the angular distribution of these photons is not described accurately by 
equation (2), even if A max /L min is very small. Indeed, the Eddington factor at energy e in is equal to unity everywhere in the slab. In 
this example, equation (2) is accurate only at energies sufficiently below the injection energy e in . 

The electron bulk velocity V is a vector quantity, and hence the relative sizes of the terms of first and higher order in V in 
the expansion of the transfer equation and its moments depend strongly on the relation between the angular dependence of 
the bulk velocity, the radiation field, and the differential scattering cross section. Therefore, even when the diffusion approx- 
imation is only slightly inaccurate (as for example when the Eddington factor differs only slightly from |), this inaccuracy 
produces additional terms in the moment equations that are of second order in V and that are therefore of the same 
magnitude as the second-order terms that would be present if the radiation field were given exactly by equation (2). 

For these reasons, we derive the radiative transfer and moment equations without making use of the diffusion approx- 
imation. 


3. ASSUMPTIONS, DEFINITIONS, AND APPROXIMATIONS 

In the sections that follow we assume that the electron gas is nondegenerate (electron occupation number <^1). For 
conciseness we shall refer to it as a fluid, without implying anything about whether it is collisional or collisionless. We assume 
that photons are scattered only by electrons, neglecting scattering by any other particles, and set h = c = fc B = 1, where h is 
Planck’s constant, c is the speed of light, and k B is Boltzmann’s constant. We indicate quantities evaluated in the rest frame of 
a particular electron by a subscript zero and quantities evaluated in the inertial frame momentarily comoving (locally) with 
the fluid, which we call the “ fluid frame,” by a subscript f The “ system frame” may be any global inertial frame (such as the 
frame at rest with respect to the center of mass of the system, if it is inertial). Quantities evaluated in the system frame have no 
subscript. 

Bulk velocity and temperature . — We define the fluid frame as the frame in which the energy and momentum flux density 
both vanish (Landau & Lifshitz 1987, p. 505). In this frame, collisions between the electrons tend to establish an isotropic 
velocity distribution. 

The system-frame three-velocity u of a given electron is related to its fluid-frame three-velocity u f by a Lorentz boost. The 
first and second moments of u in the system frame are 

<«> = V (3) 

and 

<u 2 > = <f 2 > + F 2 . (4) 

Here V is the three-velocity of the fluid as measured in the system frame, the angle brackets denote the average over the 
velocity distribution of the electrons in the fluid, and 


v = u — V 


( 5 ) 
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<u}>f = ^ + 

m 0 


<$■ 


( 6 ) 


and the second moment of the electron peculiar velocity distribution in the system frame is 

<v 2 )=^ + 0 (^V 2 ). 

tn e \m e ) (V 

ra^atioif field interacting^Rh Amoving 1 efectron^uid^We shalD-efe 6 fTtT ^ eqUation that describes the evolution of a 
written in terms of the photon mode occ^a on 1,1 t ^ aS the photon kinetic e( * uation when it is 

the specific intensity of the radiation field. The two descrintionfa ° nu iatlle transfer equation when it is written in terms of 
transfer equation is more ™ <**’ ** Miha,aS < 978 ’ * ^ but the radiative 

radiation field by the number „(/, a, of photons with energy . 
consider only unpolarized radiation). The first few moments of nffe) are^ C ependence on P olanzat 'on state, because we 


n = in \ n(l ' e ' dQ • 

n ‘ = ^ J «(/, e^dO , 

n i} s ~ J nil, e)fl j dn , 
n lJk = ~ J »(L E)Pl J fidQ . 


(8) 

(9) 

( 10 ) 
( 11 ) 

P 0 ^ 0 " a " d P^on **** been *>r brevity, 

its zeroth moment. de P™dence of the photon occupation number on / and e in order to distinguish it clearly from 

•> ^ *><« fi); here 


H 


K 


■if* 
■if 


e)dSl = 2 e 3 n , 


/(/, s)l‘dQ = 2e V , 


( 12 ) 

(13) 

(14) 

(15) 

where again we have suppressed the dependence of the moments on position and photon enerev 

(see Be' rest e tsk i E L i fs hi t z,& * P u a e vs k fi T 9 71 ! ^297^°"’ ^ differentia ‘ Cross section for scattering of unpolarized radiation is 


= in\ /f/ " £)l‘l J dn = 2e 3 n iJ , 
Q iJk = ~ [ /(/, fi)/W k dSl = 2e 3 n iJk , 


da 

dQ n 


a _ 3ff T AoV 
\ 1 671 \e 0 J 

where o, is the Thomson cross section, and I, are , he energy and direction of the incident photon, and 


7 + r- 1 +('«•«' . 


£n = 


1 - (fio/Wj 1 - /„ • /y 

and f' 0 are the energy and direction of the scattered photon. 


( 16 ) 


( 17 ) 
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The energy and d, rechon of propagation of a photon in the electron rest frame are related to the same quantities in the 
system frame by the Lorentz transformations 


eo=F (l-/-«)- ?(1 + /0 . H) 

(18) 

and 



(19) 

The photon phase-space volume is Lorentz invariant (see Mihalas 1978, p. 495), i.e.. 

(20) 

£ ch d£l = £o ds 0 dn 0 -> 

s0 

n{ly c) = n 0 (io, e 0 ) ■ 

(21) 


sySatltX^^ 

of order (efmju or higher. These equations are therefore valid when 

( 22 ) 


« 1 


and 


e 

m 0 




These conditions are usually satished in accretion onto white dwarfs and neutron stars bn, are not satisfied in accretion onto 
b ' radiahve'transfer equahons in ,h. fluid frame are correct to all orders in V bu, 
only to first order in TJm e and 8,/m,. These equations are therefore valid when 

(23) 


m e 


and 


nK \ in 


e f nr 


f<© 




and can be used where V — 1- 

4. PHOTON KINETIC AND RADIATIVE TRANSFER EQUATIONS AND THEIR ZEROTH AND FIRST MOMENTS 

4.1. Photon Kinetic Equation 


The photon kinetic equation in the system frame is 


a n(k) 




Wik'p, kp')n(k')L\ + n(k)lfip)L 1 -/(/>')] - W(kp, k'p')n{k )[ I + n(k')lfiP)U HP')Vi 


(24) 


where k" = (e, ef), n(k) is the photon occupation number ^ is the ‘wHtin^thc totSl derivative on the left 

side'we have Sh^instein ^Ifio^onvention The left side of equatio^ (24) is manifestly covanant. The right side ,s 

then much stmpler. We will therefore assume for the i^ «n ^evaSS in the system frame and the right side in the 

rS:^n f » tfle — - * Ofthe photon ,n the system frame. The resulting 

equation for electrons moving with velocity u - pj (ym,) is (cf. Peeble . P 


(<?, + /T,)en((, e) = e 0 n, 0 dS2 


, e o 


Pt-0 

d«o 


do 

dQ' 


n 0 (Joi Ep)[ 1 + n o0o’ Eo)] (^dQ'J 0 


l ^ a \ n 0 {l u , £ 0 )[ 1 + n o(("o’ 


Eo)] 


(25) 


• i .• ic tVif. flfrtron density in the electron rest frame, f 0 and 6 0 are 

where P = P et. P f = P/Px\ x‘ are the spatial coordinates, 'Vo • , , m . -rn. e r ac * or preceding the cross section in 

th^f^nerm^f^he^iSfision^ntegral 0 ^ for Se different phase spaces of the incident and scattered 

e/m,. This gives (see Peebles 1971, p. 204) 

* f . . ,r . _w r h.nj.i(l - /„ • l\\.P. nil'. s) + l r(/o - >o) • «] 2 e 2 ^ 2 n(K e) , (26) 
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where 8 e = 8/de and d 2 = 8 2 /ce 2 . We then use this result, the Lorentz invariance of the photon distribution, and (da/d£l') 0 and 
n e0 e 0 , expanded to second order in u, to obtain the approximate photon kinetic equation that describes the effects of 
scattering by electrons with velocity u. 


( 8 , + 1 ‘dMl, e) = 


3n e a 1 

I6n 


dQ'ol 1 + (4 ■ 4) 2 ] 1 - 2 — (1 - f 0 • /o) j[l - (/; • u) + (/ 0 • u) 2 - u 2 ] 
L m e J l 


+ [(4 - 4) • « + 2(4 • «X4 - 4) • «K + 2 C(4 - 4) • «] 2 e 2 ^ 2 + ^ (1 - 4 * 4X4 + ec £ 


+ 2 — (1 - / 0 • /«)«(/, «X2 + e(T) 


|n(/\ s) — n e <r x ^l — 2 “^(1 — f • «)n(f, e) • 


(27) 


In Appendix A we give the photon kinetic and radiative transfer equations that are obtained by averaging equation (27) 
over a drifting, relativistic Maxwellian electron velocity distribution (see eqs. [3]-[7]). The moment equations derived in the 
next two subsections can be obtained by computing the zeroth and first moments of the equations given in Appendix A. Here 
we follow the simpler approach of first computing the moments of equation (27) and then averaging them over the electron 
velocity distribution. 

4.2. Zeroth Moment and Radiation Energy Density 

We compute the zeroth moment of the photon kinetic equation by first integrating both sides of equation (27) over all 
directions. Making use of the Jacobian 


d(a 

d(Q 


l Qo) /«o V V \ _ fl 

! 0 , Q') \e de J\e' 0 de' 0 J ~~ 


2(/ 0 * tt) + 2(1' • u) - 2 u 2 + 3(/ 0 - u) 2 4- 3(f ■ u) 2 - 4 </ 0 • u)(f • i/) 


(28) 


we find 


) ' n + s ‘ n ‘-^ 


2 n e f da [ <K) 0 [1 + (to • /o) 2 ][ 1 - — (1 - t 0 ‘ 4)1 

n J J L m e J 


x H [* - 3/< 2 + 3(f • u ) 2 + 2 (f ■ «)] + (4 • «)[1 + 2(7' • «)]ec? £ + \ (4 • u) 2 e?? 2 + — (4 + 




- (4 • «){[3 + 6 (f •«)] + [!+ 2(f • «)M + (4 • u)E 2 dl} + (/" 


’o ■ "r( f ’ + 


+ 2 rt « 


— (l 0 • /q)( 4 4- £<3 £ ) -f (1 — l * /')«(/, e)(2 + fid £ ) 


e) - n e a T ^ 


Next we integrate over dil 0 and then over dQ\ using definitions (8H1 1). The result is 


1 


n e o T 


(d t n 4- di n l ) = (3 4- edJn'Uj 4- 


n + ^ 


36 34 


1 - 2 — |n - n l Ui 


11 


(29) 


— (4 + <*<?«) + V ( 4fi ^ + fi2 ^r 2 ) I" M ^ + ]q ed * + ™ £2(1 * * nUUi Uj ~ * 


20 


+ - — )(2n 2 - 2 nW 4- 2 n iJ n iJ - 2 n iJk n ijk + ned r n - n l ed r ri 4- n^ed r n ij - n ijh e8 t n { ^ 

2 \mj 


(30) 


where we have again used the Einstein summation convention. Finally, after averaging equation (30) over the electron velocity 
distribution, we obtain the zeroth moment of the kinetic equation that describes the effects of scattering by a fluid of electrons, 
namely, 

— (d t n 4- d f «‘) = (3 4- ed^Vi 4- — (4 4- ed £ )n 4- ^ «t ,2 > + F 2 )(4 ed E + c 2 d 2 )n 
n e <j T m e 3 


36 34 11 _ 

lo + To + 20 ' a ' 


iff - 5 •<»*>) + (1 + 15 < + 55 ^ V, - 3 " 




3 / £ 

+ 2 \ni 


1(2 n 2 — 2 n l n l + 2n lj n lj — 2 n ljk n ljk 4- ned € n — n l ed t n l + n lj ed £ n lj — n ,jk ed e n ljk ) 


(31) 


This equation is valid for any arbitrary (possibly anisotropic) distribution of electron velocities in both the diffusion and 
free-streaming regimes. 

If the electron velocity distribution in the fluid frame is a relativistic Maxwellian with temperature T e , then (see eqs. [3]-[5] 
and [7]) 

<ViVj) < L ’ 2 )^ ~ — S u , 

3 m„ 


( 32 ) 



888 


PSALTIS & LAMB 


Vol. 488 


and the zeroth moment of the photon kinetic equation can be written, to the same accuracy as that of equation (3 1 ), as 

1 


n e 


(d, n + <?; n') = (3 + ed, 


?>‘K + 


— (4 + ed c ) + 
m„ \m 


l + y ) 4e ‘ i - + + 05 + 15 “*• + S V ‘ ~ 5 " y ‘) 


3 t 


+ - ( — )(2n 2 — 2n l n l + 2 n lj n lj — 2 n ljk n ljk + nsd E n — n l ed c n l -f n lj ed c n tJ — n ljk ed E n ljk ) , (33) 

2 Vm, 


where we have used the relation = n . The corresponding zeroth moment of the radiative transfer equation, with emission 

and absorption included, is (see Appendix B) 


f + h - H( 


£ - 4T 


T V 2 

+ — e8 2 (eJ) + ed t H‘V, + — [-4 e 8 c J + eS 2 (eJ)^ 
m. 3 


+ (i^ + M - 1 •'"'0 

+ 5 - J )?~ (•»•*■ - H ) 7 - + (*«* - K “) 7 - {•*<** - e") 

H 


+ ( 1 - - Vhd c + - vh 2 dl )r, c 


{( J + \ jy2 ~ H ‘ v 7 + KiJVi Vj + \ Jy2 ~ H>Vi £dc + \ K ‘ iy ‘ ’ 


(34) 


where rj E and are the emission and absorption coefficients, which are defined in the fluid frame but evaluated at the energy e 
of the photons in the system frame. 

4.3. First Moment and Radiation Flux 

We compute the first moment of the photon kinetic equation by multiplying both sides of equation (27) by l and then 
integrating over all directions. In performing the integration we use the transformation (19) and the Jacobian (28), 
integrating first over dQ 0 and then over dQ f . The result is 

1 ( d t n l + dj n ij ) = — n* - ~ Jj- (1 + filin' + j ( — nu i + 3 n lj Uj) - ~ (8 n lJk UjU k - 2n l u j u J ) - ^ [3 ed c (nu<) + ed c (n lJ Uj)] 


n e c7j 


10 


- + 9 ed E (n ijk UjU k ) - ed r n i u 2 '] - [3 e 2 d 2 (n j u j )u i + e 2 d 2 {n ijk u j u k )'] 

+ - ( — )(2/tn i — 2n J n lj + 2 n jk n ljk — 2n Jkl n ,Jkl -f n l ed t n — n lj ed e n j + n ljk ed € n jk — n ljkl ed E n jkl ) . (35) 

2 \mj 


Averaging equation (35) over the electron velocity distribution, we finally obtain 

2 £ 

5 rn 


1 (d , ri + djti ij ) = -n‘ - ^ (1 + sd t )n‘ + | (n i} V } - 1 nV]j - ^ ed t (n iJ Vj + 3nf)) 


- j (4 n iJk Vj V k - n‘Vj Vj) - T e dJ9n ijk Vj V k + 10 n j Vj V, - n% Vj) 

- T E 2 d 2 (n iik Vj V k + 3n j Vj V, ) - 1 {4n iJ \ Vj v k ) - rt<.VjVj» 

- T ed t (n iJk (VjV k y + \0n\vjV;) - n‘(VjVj » - ^ e 2 d 2 {n ijk <VjV k ) + 3n i (v j v i )) 


3 / e 


-f- ( — h2 nn* — 2 n j n lj -f 2n jk n ljk — 2 n jkl n tjkl -f n l ed c n — n lj ed c n j + n ljk ed € n jk — n ljkl ed r n jkl ) . (36) 

2 \mj 

If the electron velocity distribution in the fluid frame is a relativistic Maxwellian with temperature T e , then 


and 


nKvjVi) = ^ «' 
m c 


n ijk <V; v k y = — ri, 

m. 


(37) 


( 38 ) 
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and the first moment of the photon kinetic equation reduces to 




1 + c - 1 + + 7 — (1 + 4*0. + 


(d, ri + d 


+ 5 r ! ‘ J ^ “ 3 nVi ) ~ TO edc(niJVj + 3nVi) ~ 5 {4n ‘ JkV J K ~ n ‘ V J V > ] 

- T Ed e (9n ijk Vj V k + 10 n j Vj V, - n% Vj) - -1 e 2 8^n ijk Vj V k + 3 n>Vj V.) 

+ 2 ~ 2 " J ” U + 2 n Jk n ijk — 2n Jk, n ijkl + n‘i:d F n — n' J ed c n J + n ijk ed c n' k - n ijkt ed c n Jkl ) . (39) 

The corresponding first moment of the radiative transfer equation with absorption and emission included is (again see 
Appendix B) 

d, //■ + dj K‘ j =Lw-~ V - - 5 ~ edf el//' + T [(9 - ed t )K% + (7 - 3 ed,)J VJ 


- 5 l ed \—) + m e Sdte f + T0 C(9 -^^ + (7 -W^ 

+ jo (_6 + 8ed * ~ V, + i (1 + sd c )H‘Vj K. - -i (l + 3ecl, + e 2 d*)Q ,Jk Vj V k 


3 / e \| H' K" O iJk RWTl 

+ 4 {-)[^ J - J ) pr - K H’ - in -p + (s<l K* - K> k ) ^3- - (ef 1 , Q» - Q*') £5- J Jn, <r T 

-[(' + 1 ^ V ) W ~ K " V > + (2 H ‘ V > V * K ' iV > + 2 QUkVj K ‘) £< - + 2 Q ‘ JkVj ^ Z. + 5 


where R lJkl is the fourth moment of the specific intensity. 

4.4. Equations in the Fluid Frame 

The photon kinetic and radiative transfer equations as well as their moments take their simplest form in the fluid frame, 
since V = 0 in this frame, by definition. These equations can be written choosing as independent variables either Eulerian 
coordinates fixed in space and time or Lagrangian coordinates comoving with the fluid. 

In terms of the fluid-frame coordinates, the transfer equation (A8) and its zeroth and first moments, equations (34) and (40), 
become, in the fluid frame, 

( d >, + ! f 3.M- E /> = »«./ ff if 1 At /t e f ) + - [se 2 J f + H‘ f + ‘J’iKj + sef Qf + <JT Rf ') 


d, f J f + d if H‘ f = 




4e 3 \m e 


+ 4? \ » - 'X-// -/>//> + , (41) 


+ t( — IT J, 


- I « r 0 f/ H‘j - H' f \ J- + ( S/ d Cf K'j - Kf -f - (e f d Cf Qf - Qf) 



where the coefficients !£ x -<£ 6 , which are given in Appendix A, are to be evaluated at V = 0. 

The left sides of the above equations can also be written in terms of the system-frame Eulerian coordinates, as in equations 
(95.9), (95.1 1), and (95.12) of Mihalas & Mihalas (1984). The resulting equations, which are often called mixed-frame equations, 
are correct to all orders in V and can therefore be used in situations in which the bulk velocity is relativistic. The velocity- 
dependent terms that appear in the mixed-frame equations arise from the Lorentz transformation of the left sides of equations 
(41H43), whereas the velocity-dependent terms that appear in the system-frame equations arise from the Lorentz transform- 
ation of the scattering integral on the right sides of these equations. 

The right side of the radiative transfer equation (41) can also be used in the formalism developed by Thorne (1981) for 
solving Comptonization problems in general relativity (cf. eqs. [6.10] and [6.13] of Thorne 1981). 
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5. DISCUSSION 

In the previous section we derived the radiative transfer equation in both the system and fluid frames, taking into account 
absorption and emission, as well as spontaneous and induced Compton scattering. In this section we first show that our 
equation reduces to the Kompaneets equation in the appropriate limits, and we call attention to several errors and misunder- 
standings in the literature. Next we discuss the moment equations for an anisotropic radiation field in a static medium and 
show that the radiation force, and hence the critical flux and luminosity, generally depend both on the photon energy 
spectrum and on the electron temperature. We then consider the equation for the zeroth moment of the specific intensity for 
moving media and show that if the radiation field is isotropic, the terms in the transfer equation that are second order in the 
electron bulk velocity produce a systematic increase in the energy of the photons that is completely analogous to the 
systematic increase in the energy of the photons produced by the electron thermal motions. We also show that Com- 
ptonization by electron bulk motion occurs whether or not the radiation field is isotropic or the bulk flow converges, and we 
give estimates for the timescales on which the photon energy distribution changes because of systematic downscattering and 
upscattering caused by the electron thermal and bulk motions. We derive a new, more general condition for determining when 
Comptonization by the electron bulk motion is more important than Comptonization by the electron thermal motions. We 
conclude by indicating how the transfer equations we have derived can be solved using the method of variable Eddington 

factors. 


5. 1 . The Kompaneets Equation 

When the radiation field is isotropic and there is no bulk velocity and no absorption or emission, equation (34) reduces to 
the Kompaneets (1957) equation, 1 


47 : 


<\.J = f- d c (eJ) - ed c J+ ^ ed?(sJ) + — (ed c J - J) 


(44) 


m 


m 


m 0 


where on the left side we have introduced the differential Compton time dt c = n.Ojdt. The first two terms on the right side of 
equation (44) describe the effects of systematic downscattering and upscattering of the photons by electrons. The third term 
describes the diffusion in energy produced by the thermal motion of the electrons. The last term describes the effect of induced 

Compton scattering. ...... , ^ 

As noted by Kompaneets (1957), any Bose-Einstein distribution is a stationary solution of equation (44). Induced Compton 
scattering cannot change the chemical potential, because it does not change the number of photons. Hence, the statement by 
Pomraning (1973, p. 193; see also Rybicki & Lightman 1979, p. 209) that the Planck spectrum (the particular Bose-Einstein 
distribution with zero chemical potential) is the only stationary solution of the Kompaneets equation if induced scattering is 
included is not correct. 

Integrating equation (44) over energy shows that if the electron temperature is equal to the Compton temperature 


T c = 


<£ 2 ) - <ne 2 ) 

VeS 


(45) 


the photon energy density remains constant although the photon spectrum may evolve with time; here the average is over 
photon energy, using the photon-number density N(e) a (J/e) as the weighting function. Note that for a given photon 
spectrum, induced Compton scattering always decreases the Compton temperature. 


5.2. Implications of the Moment Equations for Static Media 

As we mentioned in § 1, the Kompaneets equation (44) is not strictly valid for astrophysical systems, since it requires the 
radiation field to be isotropic, and hence that no radiation leave the system. If the radiation field is anisotropic, and there is no 
bulk motion, the system of equations 


d,J + (%H‘ = r, - X J + 


" fe - 4 T e \ T„ , I 
n e o T £<y 1 H £<> 2 £ j 

L V m e ) J 


(46) 


and 


. 2 [/T e -3e\ ,/e-4TA T e 1 

,, H> + ? , K" - — „ x W — - - «, J + ^ “ • ‘f 


(47) 


must be solved simultaneously together with closure relations for K li and J, which are usually introduced as variable 
Eddington factors (see Mihalas 1978, p. 157). For simplicity we have neglected the effects of induced Compton scattering in 
equations (46) and (47). The term (2n e o T /5)l(T e - 3e)/roJH' on the right side of equation (47) does not appear in the first 
moment of the transfer equation derived by Chan & Jones (1975), Payne (1980), and Madej (1989, 1991), because these 
authors did not use the appropriate relativistic scattering cross section. All the terms in the square brackets on the right side of 
equation (47) were neglected by Blandford & Payne (1981a), but this is generally unsafe. 


• In their derivations of the Kompaneets equation, Rybicki & Lightman (1979, p. 213) and Katz (1987, p. 100) did not take the different phase spaces of the 
incident and scattered photons into account in the collision integral. However, following Kompaneets, they evaluated the integral by using photon 
conservation and the thermodynamic equilibrium photon distribution, rather than by performing the integration directly, and they thereby obtained the 
correct result despite this error. These authors also used the Thomson approximation to the Klein-Nishina cross section. In general this introduces an error 
of the same size as the systematic downscattering term (see below), but this error vanishes if the photon distribution is isotropic. 
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In a static medium and in the absence of absorption, emission, and induced scattering, F, the radiation force per unit 
volume on the electrons (see also Lamb & Miller 1995) can be obtained by integrating equation (47) over photon energy, 
which gives 


F = 


4n 

c 



8 <^ + 2 
5 m e 



(48) 


where <£>„, is the average photon energy, using H‘ as the weighting function. Blandford & Payne (198 la) neglected the terms 
in the square brackets on the right side of equation (47), and, therefore, the radiation force given by their equation is incorrect 
to first order in e!m e and TJm e . Equation (48) shows that the volume radiation force on the electrons produced by Compton 
scattering of the photons is different from that obtained in the Thomson limit and that it depends both on the photon 
spectrum and on the electron temperature (see also Fukue et al. 1985). Hence, to first order in e./m e and TJm e , the energy- 
integrated critical radial radiation flux ,F cril that produces a radially outward radiation force that exactly balances the inward 
gravitational force of a massive object depends on the photon spectrum and on the electron temperature. For completely 
ionized hydrogen gas, the critical radiation flux at radius r is given implicitly by the equation 


erit 


crn p GM / _ 8 + 2 7;\ 1 

r 2 t 7 t \ 5 m e mj 


(49) 


where m p is the proton mass, M is the mass of the object, and (e)^ £Cr h is the average photon energy, using & c c Tlt as the 
weighting function. For example, in the spectral formation region of many neutron star low-mass X-ray binaries, <e)^ £Cri t is 
~ 1 ke V, and T e can be -25 keV (Miller & Lamb 1992; Psaltis et al. 1995), in which case the critical flux is - 10% smaller 
than the usual Eddington critical flux computed by assuming Thomson scattering. Note also that to this order, the critical 
luminosity L cril = 4 nr 2 ^ crit generally depends on radius, because both <(£)^ fCr ii and ^generally depend on radius. 


5.3. Implications of the Moment Equations for Moving Media 

Consider now a medium in which the electron bulk velocity is not zero. Suppose first, for simplicity, that the radiation field 
is isotropic in the system frame. To make this situation concrete, consider a thought experiment in which electrons are moving 
with uniform and constant bulk velocity V through a box with sides of length L. The electrons are assumed to be able to pass 
through the walls of the box, whereas photons are confined inside the box but have a mean free path to scattering by electrons 
that is much larger than L. Under these conditions, scattering of photons by the walls of the box, which is much more frequent 
than scattering of photons by electrons, keeps the photon distribution nearly isotropic. 

Tn this situation the equation for the zeroth moment of the specific intensity reduces to 




— djiej) + 
m e 



(- 4 ed c + ed 2 e)J , 


(50) 


where, for simplicity, we have neglected absorption, emission, and induced scattering. Equation (50) shows that when the 
radiation field is isotropic in the system frame, Comptonization by the bulk motion of the electrons is described entirely by 
terms that are second order in V ; all terms that are first order in V vanish identically. 

Suppose now that (i) photons with energies e T e + \m e V 2 are injected into the box with an isotropic momentum 
distribution, (ii) the photons are allowed to remain in the box for a distribution of residence times that decreases exponentially 
with the residence time, and (iii) 



(51) 


where f res is the mean residence time measured in units of the Compton time (\jn e <r r ). The spectrum is then a power-law with a 
high-energy cutoff (seeSunyaev & Titarchuk 1980). At high and low energies, the spectrum is 


where 


J = 


e exp [-6/(7; + ±K 2 )] , 


£<T e + X V 2 , 
e > T e + jni,, V 2 , 


(52) 


a = 



(53) 


and the factor e 3 arises because we are considering the energy density rather than the photon occupation number. This 
solution is a simple generalization of the solution to the Kompaneets equation obtained by Shapiro, Lightman, & Eardley 
(1976; see also Rybicki & Lightman 1979, chap. 7). It shows that the terms in the transfer equation that are second order in the 
electron bulk velocity produce a systematic increase in the energy of the photons that is completely analogous to the systematic 
increase in the energy of the photons produced by the electron thermal motions, as argued in § 2.2. 

When the photon momentum distribution is not perfectly isotropic, the characteristic timescales on which the photon 
distribution changes because of systematic downscattering and Comptonization by the electron thermal and bulk motions are 
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(see eq. [34]) 


down ~ < 

(54) 

( 4T <\ 

h h ~ « e ff T 1 , 

(55) 

\mj 


ty' ~ n * ff T y 'yjJ ’ 

(56) 

f 4V2 \ 

t v z ~ y\ e ^ J ’ 

(57) 


where the last two timescales have been estimated from the terms that are first and second order in V. 

Comparison of rates (55) and (57) shows that Comptonization by the electron bulk motion is more important than 
Comptonization by the electron thermal motion if 

V 2 > IZf (58) 

m e 

In fact, as explained in § 2.2, Comptonization by the bulk motion occurs whether or not the radiation field is isotropic or the 
bulk motion converges (see also eq. [50]), contrary to the impression given by Blandford & Payne (1981a). Comparison of 
rates (56) and (57) shows that ifJ is sufficiently large compared to H\ the effects of bulk Comptonization described by the terms 
that are second order in V are dominant compared to the effects described by the terms that are first order in V. 

In estimating the characteristic timescale for bulk Comptonization, Blandford & Payne(1981a; see also Blandford & Payne 
1981b and Payne & Blandford 1981) used only one of the several terms in their equation that are first order in V, neglecting 
other terms that generally also produce systematic upscattering or downscattering of photons. As a result, the characteristic 
timescale that they obtained for bulk Comptonization is proportional to (V • F)/3, which they assumed to be proportional to 
n e a T V 2 / 3. This assumption is not generally valid. When it is, their expression suggests that bulk Comptonization is less 
important than thermal Comptonization if V 2 < 12 TJm e . However, comparison of rates (55)-(57) shows that the bulk 
Comptonization terms that are second order in V are already as important as the thermal Comptonization terms when 
V 2 ~ 5TJm e , and hence, that these terms can be more important than the thermal terms even if V 2 < \2TJm e . 

5.4. A Method of Solving the T ransfer Equation for Moving Media 

In more realistic models, the radiation field is not isotropic in the system frame in the presence of bulk motion (see, e.g.. 
Miller & Lamb 1993, 1996), because electron scattering tends to isotropize the photon distribution in the fluid frame. In this 
case, the full radiative transfer equation must be solved. This can be done by using the method of variable Eddington factors 
(Mihalas 1980; see also Mihalas 1978, p. 201), in which the radiative transfer equation and its zeroth and first moments are 
solved iteratively. In this approach, the second and higher moments of the specific intensity are related to the zeroth and first 
moments via variable Eddington factors. The zeroth and first moments of the radiative transfer equation are then solved by 
using initial guesses for the variable Eddington factors, and the source function is computed from the calculated moments of 
the specific intensity. The radiative transfer equation is then solved, the Eddington factors are updated, and the whole 
procedure is repeated until convergence is achieved. 

A detailed study of the solutions of the equations derived here for realistic models of astrophysical systems will be reported 
elsewhere. 

It is a pleasure to thank G. Baym, V. Kalogera, D. Mihalas, J. Poutanen, L. Titarchuk, R. Turolla, L. Zampieri, and S. Zane 
for helpful discussions and comments on the original manuscript. We are also grateful to V. Kalogera for help in checking the 
derivations. This research was supported in part by NSF grants AST 93-15133, AST 96-18524, NASA grant NAG 5-2925, and 
several Rossi X-Ray Timing Explorer grants. 

APPENDIX A 

PHOTON KINETIC AND RADIATIVE TRANSFER EQUATIONS FOR SCATTERING BY AN 

ELECTRON FLUID 

In this appendix we give the photon kinetic equation and the corresponding radiative transfer equation for scattering by an 
electron gas with temperature T e and bulk velocity V. 

We start from the photon kinetic equation (27) for scattering by electrons with a given velocity, transform all the quantities 
to the system frame using equations (18) (19), and then average over the electron velocity distribution. The resulting photon 
kinetic equation is 

-i- (d, + /Y~,)n(f, e) = i# A n(l, s) + l (# 2 n + + Mitt* + 3tfn ,ik + #f l ri Jkl ) 

", 4 

+ - ( — ]«(/, eK2 + edjn - Iff + l'lff j - l'Pl k n ijk ) , ( A 1 ) 

2 \mj 
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where 


m x = -1 +2 — + (/• V) H 
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01 2 = 1 + 2 — + — ed c + 2 — (/ • V ) + (/ ■ V ) 2 - V 2 + e 8 £ 4 - (/ • V ) + (/ • V ) 






L m e 


2 ] + \ 


2 — + (/ • F) 3 


= -2 — /' - — Ed , l ‘ - 4 — 1 ‘ + 2 Vi - 2(1 • F)/ f - 2(1 • F) 2 / 1 + 2 V 2 1 ‘ 
nr , m „ 


+ sd e -4 T f + y . _ 4(/ • F)Vj + 2(f • F) 2 / 1 

W 1 „ 


I £ d? 
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2 — I 1 - 2(1 • V )Vi 
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']■ 


m = l ‘ l j + 2 — I'V + — ed c Pl j - 6 T i‘[J + d . K)/7 j - 2/‘F, + (/ • VflH 1 - 8 (Z • F)/'F - 3 V 2 l ‘ l J + 4K, F 


+ 4 — W - (l ■ V ) l ‘ l J + 3 V i V ,-( t - V ) 2 il } 
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e5, W + 4 — /'7 J 7* + 4/7'F* - 8/T 7 F t + 4 (f • F)/7 J 'F* 


2 — ri j + v i v i + ( i - vfiv 

m„ 


"7 J , 


'>;■ K»1 + 1 6 2 0 2 |^ 


- 2 — Z‘7 J 7* — 2(/ • F)WF* 

m e 


"/'KJ , 


= -2 — /7'/ k - — ed, /7 J 7* + 4 — 

m t , m e m 

+ fid -4 T /7 J /‘ + /7 J K - 4(/ • F)/7 y K - 21% V k 
L m e 

' = lOl‘l*V k V, + 5ed e lW k v i + \ e 2 ^ 2 /‘T'X F,. 

The corresponding radiative transfer equation is 

— (5, + /*&)/(£ e) = if, /(/, s) + 7 (,r 2 ./ + S/ y \ II‘ + . r S'jK ,J + 2 ffQ ijk + * ¥f l R iikl ) 
n e o t 4 

+T ( — )/(/“ fiXfi^v - w - i‘h' + ri J K ij - m k Q ijk ) , 

4£ 


where 


= -1 + 2 — + (/-n, 
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^2*' = l‘l% V, + 2fifi>, l i l i V k V, + - c 2 cf l‘l J V k V, . 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A 13) 


(A14) 
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EMISSION AND ABSORPTION PROCESSES IN THE SYSTEM FRAME 


In the absence of scattering, the transfer equation in the system frame becomes simply (see, e.g., Mihalas & Mihalas 1984, p. 
422) 


(d, + l%m e) = 






(Bl) 


where again the subscript /refers to quantities in the fluid frame, and we have, for simplicity, assumed that the emission 
coefficient r\ and the absorption coefficient x are isotropic in this frame. In writing equation (Bl) we have also used the fact 
that the absorption and emission coefficients are defined in the fluid frame and that (rj/e 2 ) and (ex) are Lorentz invariants. 

We can now expand rj(e f ) to second order in V as 

nU'j) = >/[£ + (£/- £)] + (a f - e)d c ri L + j(a f - e) 2 3 2 ri c ~rj e + ({ V 2 - 1 ■ V)ed c ri c + (/ • V) 2 e 2 d 2 n , , (B2) 

where rj t . = rj(e). We also expand x(£f) in a similar way. The transfer equation to second order in V then becomes 

(/>, + /V',)/(/, £) = {[1 + 2 (/ • V) + 3(/ • V) 2 - V 2 \ + [-(/• V) - 2(7 • Vf + |K 2 M + W ‘ V) 2 e 2 tfK. 

- {[1 - it’ V) + j V 2 ] + [-(f • V) + (/ ■ V) 2 + + i(/- V)hi \} X J(l £) • (B3) 


We can then obtain the zeroth and first moments of the transfer equation. 




- iJ + + £ c \ x £ ] - lv 2 e S t Ti c + iV 2 E 2 d 2 ri c 

\x, -iv 1 -{{ -IV 2 + K iJ V ( Vj)Ed, :Xc - \K' J V t V jE 2 d 2 Xc 


(B4) 


and 


r t H l T PjK" = -Xr W + \v\2 If. - ed'fit) + + sd.xJ ~ h*U'V 2 - (\H l V 2 + h fc£ 2 ^ 2 


(B5) 


Blandford & Payne (1981a) added a photon source term to the right side of their transfer equation without adding any 
corresponding absorption term. This is fundamentally inconsistent with thermodynamics: if absorption is not included, the 
photon source can never come into thermodynamic equilibrium with the radiation field. It is equivalent to including the terms 
on the right sides of equations ( B3)— ( B5) that involve the emission coefficient but neglecting all the terms that involve the 
absorption coefficient x?.- 

Neglecting the absorption terms but not the emission term is a valid approximation only if the specific intensity of the 
radiation field is negligible compared to the source function f/\ = t\JXr. at a ^ photon energies , i.e., only if self-absorption is 
never important at any energy, anywhere in the system. This is rarely the case in astrophysical systems. For example, it is not 
the case in the Comptonizing regions around accreting neutron star and black hole X-ray sources. 

Neglecting the absorption terms in the transfer equation leads to equations that have a different mathematical character 
from equations (B3) (B5), because in this case the right sides of the transfer and moment equations do not depend explicitly on 
the radiation field. 
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